library(bigMap)
source('~/bigMap/bdm_xtrs.R')
bdm.mybdm('~/omics/10xmouse/umap')
## [1] "~/omics/10xmouse/umap/"

UMAP. Outputs

omics.class.legend(ncol = 8)

load('~/omics/10xmouse/umap/10x1M_nn50.RData')
load('~/omics/10xmouse/umap/10x1M_nn100.RData')
load('~/omics/10xmouse/umap/10x1M_nn200.RData')
load('~/omics/10xmouse/umap/10x1M_nn400.RData')
load('~/omics/10xmouse/umap/10x1M_nn800.RData')
load('~/omics/10xmouse/umap/10x1M_nn1600.RData')

Precision (beta) qMaps

k-ary neighbourhood preservation

omics.kn.plot(list(m.50, m.100, m.200, m.400, m.800, m.1600))

Running times

m.list <- list(m.50, m.100, m.200, m.400, m.800, m.1600)
r.times <- sapply(m.list, function(m) round(unlist(m$t) /60, 0))
r.times <- rbind(r.times, apply(r.times, 2, sum))
rownames(r.times) <- c('betas', 'ptsne', 'total')
colnames(r.times) <- sapply(m.list, function(m) m$Xbeta$ppx)
knitr::kable(r.times, caption = 'Computation times (min)') %>%
  kable_styling(full_width = F)
Computation times (min)
50 100 200 400 800 1600
betas 0 0 0 0 0 0
ptsne 5 6 7 11 18 38
total 5 6 7 11 18 38

Discussion